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Abstract 

We study the nuclear medium effects and nonisoscalarity correction in the extraction of weak 
mixing angle siv?6w using Paschos-Wolfenstein(PW) relation. The calculations are performed 
for the iron nucleus. The results are discussed along with the experimental result inferred by 
NuTeV collaboration. The nuclear medium effects like Fermi motion, binding, shadowing and an- 
tishadowing corrections and pion and rho meson cloud contributions have been taken into account. 
Calculations have been performed in the local density approximation using a relativistic nuclear 
spectral function which includes nucleon correlations. These studies may be useful for the ongoing 
MINERVA experiment as well as for the proposed NuSOnG experiment. 
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I. INTRODUCTION 



MINERVA [1] is presently taking data using neutrinos from NuMI Lab., and their aim is 
to perform cross section measurements in the neutrino energy region of 1-20 GeV and with 
different nuclear targets like helium, carbon, oxygen, iron and lead. Among the various goals 
of the MINERz/A experiment, one of the goal is to measure d ^q5. in the deep-inelastic scat- 
tering (DIS) region using these nuclear targets, obtaining the ratio of the structure functions 
between the different target materials and also to study structure functions in the DIS and 
transition region {2, 3|. Neutrino Scattering On Glass (NuSOnG) jj] is another experiment 
proposed at Fermi lab to study the neutrino and antineutrino charged current deep-inelastic 
scattering events to precisely measure the structure functions F% , xF^, etc. Their plan 
is also to measure sin 2 9w from z^-nucleon scattering using the Paschos-Wolfenstein(PW) 
relation j^j, a similar type of analysis that was performed by the NuTeV group js-8]. 
Recently, we have studied nuclear medium effects on the electromagnetic structure func- 



tion Ff M (x, 



l2 ) (Of, and the weak structure functions F 2 (x,Q 2 ) and F s (x,Q 2 ) 0, For 



pEM 



x,Q ) [9(, we 
results from JLab 



bund that our results are reasonably in good agreement with the recent 
121 ] as well as with some of the older experiments like SLAC [13). In 



the case of v{y) deep inelastic scattering induced processes, the results were compared with 



the available data from NuTeV, CDHSW and CHORUS experiments |14Hl6| for the weak 
structure functions F 2 (x,Q 2 ) and xF 3 (x,Q 2 ) in iron [10| and lead [ll|, as well as with the 
results of the differential scattering cross section In this work, we study the effect 

of nuclear medium and nonisoscalarity correction in extracting sin 2 9w using PW relation 
taking iron as the nuclear target. 



Paschos and Wolfenstein 



5| suggested that the ratio of neutral current to charged current 



neutrino interaction cross sections on nucleon targets may be used to measure sin 2 6w- 



R 



PW 



a{v„ N Up X) - N ^ X) _ 1 



sin 2 fl 



w 



(1) 



cr(z^ N -> fi- X) - a(u^ N -> fi+ X) 2 
The above relation is valid when there is no contribution from heavy quarks, neglecting the 

charm quark mass, assuming the strange quark and anti-strange quark symmetry, consider- 
ing no medium effect and no contribution from outside the standard model. This equation 
is valid for both the total cross section a as well as differential cross sections -r-§-, because 

ax ay ' 

the neutral current (NC) differential cross section can be expressed in terms of the charged 
current (CC) ones, and these cancel out in the quotient [l7|. For the total cross sections, 



the Paschos-Wolfenstein (PW) relationship is true under more general assumptions (see for 
instance 0, [l7] and 18]). 

The above relation is also valid for an isoscalar nuclear target for both the total cross 
sections and differential cross sections and the above Eq. ([1]) may be written as 

<r(i/„ A -> X) - tripp A-tVpX) 



R 



PW 



o(y^ A -> fi- X) - a{v^ A -)• X) 



(2) 



where cx(i^(i^) A 



VnW, 



' M ) X) is the neutral current induced neutrino (antineutrino) cross 



section, aiy^ip^) A — > u~(u + ) X) is the charged current induced neutrino (antineutrino) 
cross section for a Z=N nuclear target A. The condition of pure isoscalarity includes the 
requirement of the cancellation of different strong interaction effects which also include the 
nuclear medium effects in the ratio of the neutral current to the charged current scattering 
cross sections. 

NuTeV Collabn. jfj-s} using an iron target has measured the ratio R of neutral current 
to charged current total cross sections, used r the ratio of charged current antineutrino to 
neutrino cross sections i.e. r 
mixing angle sin 2 6 W using Eq. (J2]) as 



4^ 4-r^= v\ as h, and obtained the value for the weak 

a(i>fi A— >fi X) 2' 



R 



PW 



a(u^ A->i^ X) 
ct(V m A-^fM~ X ) 



afuft A-^v^ X) cr(u^ A-+n+ X) 
ct(P m A^fi+ X) A-tp- X) 



1 - 



R v - r R? 



where R u exp 



0.3916 ± 0.0007 and R D exp 



of sm 2 6 w is 0.2277 ± 0.000- 



o-(P M A^n+ X) 

o-(z/ M A^fJL- X) 



0.4050 ± 0.0016 



a a. 



(3) 



The reported value 



which is 3 standard deviations above the global fit of 



sin By/ = 0.2227 ± 0.0004 [19( and this is known as NuTeV anomaly. To resolve this 
anomaly, explanations within and outside the standard model of electroweak interactions 
have been looked for 20]- 36| . 

Paschos and Wolfenstein j^J relation is valid for an isoscalar target while iron is a non- 
isoscalar target(N=30,Z=26), therefore, nonisoscalar corrections are required. Furthermore, 
nuclear dynamics may also play an important role in the case of neutrino nucleus scattering. 
Various corrections made by the NuTeV Collaboration have been discussed in literature, 

n r 

but still the reported deviation could not be accounted for [8]. Theoretically, Kulagin [25] 
has investigated the effect of nuclear medium on the PW ratio and pointed out that the 
shadowing effect being a low x and low Q 2 phenomenon is small in the Q 2 region of NuTeV 



experiment [37] and observed the effects of Fermi motion and binding energy correction also 



to be small but isoscalar correction to be significant. Kumano [21| in a phenomenological 
analysis pointed out the difference between the nuclear effects in the valence u and d quark 
distributions may be a reason for this anomaly, however, this effect is too small to explain 
the anomaly. Recently Thomas |36j] has discussed various possible corrections and con- 
cluded that charge symmetry violation and the isovector EMC effect together may explain 
this anomaly. In the present work, we have used the results of our study of nuclear medium 
and nonisoscalarity correction on the weak structure functions and the differential scattering 



cross sections 
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111 ], to study these effects on the extraction of sin 2 9w using PW relation. 



We have obtained the modified PW relation for a nonisoscalar target. This study has been 



performed using a relativistic nucleon spectral function [38|, |39], which is used to describe 
the momentum distribution of nucleons in the nucleus and we define everything within a 
field-theoretical approach where nucleon propagators are written in terms of this spectral 
function. The spectral function has been calculated using the Lehmann's representation for 
the relativistic nucleon propagator and nuclear many body theory is used for calculating 
it for an interacting Fermi sea of nuclear matter. The local density approximation is then 



applied to translate these results to finite nuclei 



40 



4=11 ] . The contributions of the 



pion and rho meson clouds are taken into account in a many-body field-theoretical approach 



which is based on Refs. 



421 ] . We have taken into account the target mass correction 



(TMC) following Ref. [43], which has a significant effect at low Q 2 , and at moderate and 
high Bjorken x. To take into account the shadowing effect, which is important at low Q 2 
and low x, and which modulates the contribution of pion and rho cloud contributions, we 
have followed the works of Kulagin and Petti 4J, |45] . All the formalism is same for neutral 
current scattering as done in the case of charged current neutrino/antineutrino induced re- 
actions. For the numerical calculations, parton distribution functions for the nucleons have 
been taken from the parametrization of CTEQ Collaboration (CTEQ6.6) 

We proceed as follows. In Sect. Ull we present the formalism where we write the expression 
for the ^(^)-nucleon differential scattering cross section in subsection- IH Al the expressions 
for z/(z/)-nucleus differential scattering cross section for the isoscalar as well as nonisoscalar 
nuclear targets are given in subsection- IH Bl In subsection- IH CI we explicitly show the con- 
struction of the nuclear hadronic tensor for non-symmetric nuclear matter and in subsection- 
III Dl the nuclear corrections to PW ratio is presented. In Sect. IHH we present and discuss 



46]. 



the results of our calculation. In Sect. HVj we conclude our findings. 



II. FORMALISM 

A. Deep Inelastic Neutrino Nucleon Scattering 

The cross section for the charged current (CC) neutrino interaction with a nucleon target 

i.e. 

vt{k) + N(p) r(k') + X(p'), 1= e, n, (4) 

is given by: 

= 1 2m v 2M f d 3 k' 2rm yr f <$p\ yr (2M[\ yr f 1 \ vv 2 
° ' v rel 2E u (k)2E(p)J (2^)3 2 ^(k')iV (2nr\j{2ElJ^ b \2u J 'J^^ 1 1 

x ^^(p + k-k'-J^p^j (5) 

where / stands for fermions and b for bosons in the final state X. The index i is split into I 
and j for fermions and bosons respectively, 

T is the invariant matrix element for the above reaction and is, written as 



After performing the phase space integration in Eq. fl5]), the double differential scattering 
cross section evaluated for a nucleon target in its rest frame is expressed as, 

" Lff |K| / ™W \ raft TJ7 iV 



, (7) 



(2tt) 2 |k| V9 2 -^w. 

where Gp is the Fermi coupling constant, my/ is the mass of the W boson, /(= e, p) is a 
lepton, k is the incoming neutrino four momentum and k! is the outgoing four momentum 
of the lepton, q = k — k' is the four momentum transfer and Q',E' refer to the outgoing 
lepton. N is a nucleon, X is a jet of n hadrons consisting of fermions(f) and bosons(b) in 
the final state labeled by / and j in the following. 

The lepton tensor for antineutrino(neutrino) scattering L al3 is given by 



5 



and the hadronic tensor W^ R is defined as 



SjV X Si 1=1 



d 3 Pi TT f 2M l^ t r / 1 



(27T) 



2w 



(9) 



i=l 



where s^v the spin of the nucleon and the spin of the fermions in X. In the case of 
antineutrino (X\J a \N) is replaced by (X\J^\N). 



The most general form of the hadronic tensor W^g is expressed as 



W. 



N 

a/3 



QaQ(3 



1 



p.q 



p.q 



T v[y) 



2M 2 



e a ^P P q a wf ] + 



(10) 



where M is the nucleon mass and are the structure functions, which depend on the 
scalars q 2 and p.q. The terms depending on W4, W5 and We in Eq. ( llOp do not contribute 
to the cross section in Eq. (j7J) in the limit of lepton mass mi — > 0. 
In terms of the Bjorken variables x and y defined as 

Q 2 



x ~2Mu' y ~ e; q2 



2 p-q 

-q , v = — 



(11) 



and expressing W t N in terms of dimensionless structure functions F^' u (x,Q 2 ) 

Ff\x 1 Q 2 ) = MWf\v,Q 2 ) 
Ff\ Xl Q 2 ) = vWf\v,Q 2 ) 
F^\x,Q 2 ) = uW^\u,Q 2 ), 



(12) 



the expression of the differential cross section, for deep inelastic scattering (DIS) of neutrino 
with a nucleon target induced by charged current reaction is given by: 



dx dy tt(1 + Q 2 /M 2 V ) 2 



G\ME V 



/, m f \ /-, Mx. 



F 2 (x,Q 2 ) (13) 



± 



xy(l ) — 



F 3 (x,Q 2 ) 



In F3, +sign(-sign) is for neutrino(antineutrino). F% and F2 are related by the Callan-Gross 
relation leading to only two independent structure functions F 2 and F 3 . For l=e,/i we take 
mi = and assume Q 2 « M^- 



The nucleon structure functions are determined in terms of parton distribution functions 
for quarks and ant i- quarks given by 



F% p = 2x[d(x) + s(x) + u(x) + c(x)] F% p = 2x[u(x) + c(x) + d(x) + s(x)] 

F™ = 2x[u{x) + s(x) + d{x) + c(x)} F~» n = 2x[d(x) + c(x) + u(x) + s{x)\ 

xF^ p = 2x[d(x) + s(x) — u(x) — c(x)] xF 3 ™ = 2x[u(x) + s(x) — d(x) — c(x)} 

xF* p = 2x[u{x) + c(x) - d(x) - s(x)} xF* n = 2x[d{x) + c(x) - u(x) - s(x)} (14) 

For the neutral current (NC) induced reaction 

ui(Pi)(k) + N(p) u^k') + X(p'), 1= e, fi,r (15) 

the expression of the cross section fTl3l) modifies by changing Mw — > Mz, the mass of Z° 
boson and the corresponding NC structure functions are given by 

F^ c (up, up) = 2x((u 2 L + u 2 R )[u + c + u + c] + (d 2 L + d 2 R )[d + s + d + s]) (16) 
xF^ c (up, up) = 2x((u 2 L -u 2 R )[u + c-u-c} + (d 2 L -d 2 R )[d+s-d- s]) 

for the proton target and 

F^ c \un,un) = 2x((u 2 L + u 2 R )[d + c + d + c] + (d 2 L + d 2 R )[u + s + u + s]) (17) 
xF^ c (un, un) = 2x((u 2 L — u 2 R )[d + c — d — c] + (d 2 L — d 2 R )[u + s — u — s]), 

for the neutron target. Here ul = \ — Isin 2 ^, u R = —^sm 2 9 w and d^ = — \ + 
\ sin 2 9w, dn = \ sin 2 6 W . 



B. Deep Inelastic Neutrino Nucleus Scattering 

When the reaction given by Eq. (jl]) takes place in a nucleus several nuclear effects have to 
be considered. One may categorize these medium effects into two parts, a kinematic effect 
which arises because the struck nucleon is not at rest but is moving with a Fermi momentum 
in the rest frame of the nucleus and the other one is a dynamic effect which arises due to 
the strong interaction of the initial nucleon in the nuclear medium. For details please see 



the discussion given in Refs. 
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4l|. 



The expression for the differential scattering cross section from a nuclear target A is 
similar to Eq. ([7]) and given by: 



da, 



v(p)A 

cc 



Gp 



k! 



dE' dtt' (2tt) 



rn 



w 



rn 



2 I • L u,u W at 3 
W 



(18) 



where is given by Eq. flH]) and W^i the nuclear hadronic tensor is given by: 



W. 



v(p)A 



QaQ/3 

— ^ ~ 9o 



W^ )A (P A ,q) 



+ 



M\ 



P 



Aa 



r 



p 



P A -q 



A3 



q,)W^ )A (P A ,q) 



u(u)A 



PA,q) 



(19) 



where Pa is the momentum of the nucleus A. 

In the local density approximation the nuclear hadronic tensor W^r can be written as a 
convolution of the nucleonic hadronic tensor with the hole spectral function. For symmetric 



nuclear matter, this would be 
W v a f )A = 4 / ,/ V 



10|: 



d 3 p M 



dp S h (p ,p,k F (r))W: { ; )N , 



(20) 



(2tt) 3 E(p) J_ c 

where kp(f) is the Fermi momentum for symmetric nuclear matter, depending on the density 
of nucleons in the nucleus, i.e. kp{f) = ( ^f-p(^) ) • Sh(p°, P, kp(r)) is the hole spectral 
function, fi is the chemical potential and have been taken from Ref. Q. W"W N is the 
hadronic tensor for the free nucleon target given by Eq. (jTUI) . M and E(p) are respectively 
the mass and energy of the nucleon. 

The natural extension of the above expression for taking into account the non-symmetric 
nature of the target nucleus would be to consider separate distributions of Fermi seas for 
protons and neutrons. The expression is 111 ]: 

d 3 p M fv P 



w; 



v{v)A 



2 / d 3 r 
2 I d 3 r 



(2tt) 3 E(p) 
d 3 p M 



dp°sr ton (p°,P,k F , p )-w c 



v{v)p 



oo 

/in 



dp° sr trcn (p\pMn) -wf )n 



(21) 



(2tt) 3 E(p) "» - "<# 

where the factor 2 in front of the integral accounts for the two degrees of freedom of the spin 
of the nucleons. These are two different spectral functions and normalized to the number of 
protons or neutrons in the nuclear target. k FjP = (37r 2 p p ) 1//3 (kp :n = (37T 2 p n ) 1 ^ 3 ) is the Fermi 



S 



momentum of proton (neutron). For the proton and the neutron densities in iron, we have 
used two-parameter Fermi density distribution and the density parameters are taken from 
the Ref. (47 1. 



C. Construction of the nuclear hadronic tensor for non-symmetric nuclear matter 

The natural invariant quantities for deep inelastic scattering (DIS) of neutrinos with 
nuclei are: 

XA = 2p- g > VA= m7 ; VA = —k (22) 

where xa is the Bjorken variable in the nucleus and xa £ [0, 1]; Da is the inelasticity. These 
two variables are related to the nucleonic ones via: 

x * = -£ VA = ^ = y (23) 

where x is defined in Eq. ( TlTT) . We can see that x G [0, A|, although for x > 1 the nuclear 
structure functions are negligible. The variable ua varies between the following limits: 

< y A < - \ lAXA « — 1 - m , (24) 

so, for sufficient high neutrino energy we have < < 1. 

If we express the differential cross section with respect to these variables (xa, Va), we obtain 
the following expression in terms of the nuclear structure functions: 



dxA dyA 7T \Q 2 + m 



y\x A Fi 



l-VA- ^|^} F^ A ± x AVA (l - V f) Ff^ (25) 

For the neutral current induced neutrino interaction, the form of the differential cross section 
is the same as for the charged current induced process but with the following changes: 

U °CC . u °NC 



dx A dy A dx A dy A 
m w — > m z 

pec f nc (26) 

First we look at the denominator of the Paschos-Wolfenstein relationship, an expression 
similar to Eq. ()2]), but in terms of the differential scattering cross section, for which we 



subtract the charged current antineutrino-nucleus cross section from the charged current 
neutrino-nucleus cross section and obtain the expression as: 

d 2 a& &<il A c _G\M a E v 



dx A dy A dx A dy / 



y\x A (F? A -F? A ) 

^ _ M^U | {F „ A _ F , A) + XAyA ^ _ V_A} {F ,A + puA 



(27) 



Since Q 2 « in the present study, we have neglected the W-boson propagator term. 

We need to relate the nuclear structure functions F A to the nucleon ones via an integral 
with the spectral function. Therefore, we introduce the following notation to avoid writing 
every time the integration symbols. For example, we may rewrite Eq. f[2"Tj) with the following 
notation: 

W v f A = ( W v f )p ) + ( W v f )n ) (28) 

<*P \ °tP I proton ^ \ a/3 / s „eutron V ' 

where ( ) stands for proton: 

^ mp J ^proton 

= 2 / " 3r J (0si^y £>° sr m (p",pM r x { ; ,p (29) 

and for the neutron the expression is the same when the indices for the proton are replaced 
by the neutron indices. 

Now, if we take the three momentum transfer q along the Z-axis, i.e, g M = (q° ,0,0, q z ), 
and writing the xx-component of the nuclear hadronic tensor (Eq. [T9l) . we get it in terms of 
the nuclear structure function F\ i.e. 

W u{D)A = W^ (V)A = Fl ( } ^ = W u(D)A (30) 

Similarly, if we take the xx-components of the nucleonic hadronic tensor given by Eq. ( jTUl) . 
but remembering the fact that nucleons in the nucleus are not at rest, the xx-component of 
the nucleonic hadronic tensor is not related only to the nucleon structure function F\, but 
it is a mixture of F\ and F2 components like the expression written below: 

w v W N = w u{p)N 4. Px w u{D)N = 1 ¥Zl + Px 2 ^ Xn > rsi 1 

xx 1 M 2 2 M M 2 v N K ' 

where N = p,n; x N = £ = 2(p o g g pzt?z) and v N = £ = 
Now using Eqns. (!30l and ( 1311) we may write: 

Ff )A {x A ) _l F^{x N ) pI F^ p (x n ) \ l Ff )n {x N ) pl Ff )n {x N ) \ 

M A \ M M p-q \ M M p-q 1 ) 



10 



The difference F{ A - F[ A that appears in Eq. (1271) may then be written as: 



M A M A \M V 1 v ' 1 v u M(p-q) 



si 



+ < ^ (Fr(x N ) - Fr( XN )) + j^-^ (f^m - FtM) > (33) 



where 



F-y F^ — d v u v 

F^ v — F^ v = 2x N (d v — u v ) 

pun _ pun = _ dv+Uv = _ (pvp _ p^ 

pun _ pun = ^ +Uv) = _ (p»P _ p^ 

Here u v and d v are the valence PDFs and we are working in the so-called up and down 
quarks approximation, where we neglect strange and charm quarks contributions. 

In the case of symmetric nuclear matter, we may relate the Fermi momentum with the 
baryon density via k F = ^-p, where p is the baryon density. For a non-symmetric nuclear 
matter, we have different densities for protons and neutrons and corresponding to those, we 
also have different Fermi momenta for protons and neutrons. These are related by 

k 3 Fjt = 37T 2 p P (r); k Fjn = 3n 2 p n (r) (34) 

Instead of discussing in terms of neutron number(N) and proton number(Z) as independent 
variables, we define two independent variables A=N+Z and their difference <5=N-Z such 
that: 

N= * + l Z =*- 5 - (35) 
2 2 2 2 y 1 

Dividing the above equations by the nuclear volume V, we obtain the densities of neutrons 
and protons: 

p S p 5 . . 

Pn = - + ^; Pv = l 2 -l N (36) 

where p = p p + p n and ~ = = p n — p p . Replacing the densities for neutrons and protons 
by their corresponding Fermi momenta one has kp tP / n in terms of k F , 5 and V i.e. 



' :; 37T 2 8 3 3 3"7T 2 5 

^f, p = ™f 2~ y'i k Fn = k F H (37) 
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A non-zero value of 5 would imply that we are looking for the deviations from the isoscalarity. 
For kp, p and k F>n , we are going to perform an expansion in powers of the parameter y, and 
retaining the first order term only with the assumption that the higher orders would be 
negligible. For instance, the expansion for the Fermi momentum of the proton would be: 

and the proton's spectral function may be written as: 

7T 2 ft fl^trfi n k) 

sr ton (P°, p, k F>p ) ~ s p h ip°, P , k F - ~ s p h ip°, P , k F ) + a " l/ • 1 ' 



2Jfc|V' fcV ^ " dk 



TV 



2 (5 



^ (39) 



With a change of 5 — > —ft one gets the neutron spectral function. 

Looking back at Eq. ([27)1 . we see that the structure functions F 1 and F 2 appear as a 
subtraction and we may write F" A — F^ A as 



1 1 \ M Mip ■ q) y 2 2 ' 



si 



' F vn — F vn r> 2 



If we look at the convolution with the proton spectral function (first term of the right hand 
side), we may write this as: 



where p, p = M+ -^f is the chemical potential, which must be expanded around the isoscalar- 
ity condition, and results as: 



u_ = M+ -^~M+ \ki (42) 

Pp 2M 2M \ kpV J K J 

Inserting Eq. (1391) in Eq. (j4"T|) . we obtain the following expression with the spectral function: 



„(-fiF)} G ^ (p '^» < 43) 



| aSg(p°,p,<:) 



where 



<?*-*(p, x w ) = ( + _g_ (fr _ ) ) (44 ) 
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Integrating the first term of Eq. (H3l) . over the variable p°: 



k 2 

M+-£- 



2M 2Mk F V _ 

dp S p h (p ,p,k F )G vp -^(p ,p,x N ) 



dp u S p h ( P \p,k F )G^(p u ,p,x N ) 



-oo 

„ e!_a 

r 'ZMkrp V 



+ / * dp°S p h (p ,p,k F )G» p - pp (p°,p,x N ) 

+ S p h (fx, p, k F ) G^ p (p° = n, p, x N ) f-^-^j (45) 

k 2 

where fi = M + ^ and we have used the following property: 

/•fi+Sx 

/ dy f(y) ~ f(fi) x Sx (46) 



for — > 0. 

Let us analyze what we have obtained from the first term of the integral in expression 
1431 We have obtained an isoscalar piece (the first term in Eq. |4"5]) . which still has to be 
integrated over p°, plus a correction proportional to y. This correction would be zero for 
an isoscalar target. 

But in expression 03], we still have another integral (the one that goes with the partial 
derivative) which is already of order O (y)- Therefore, we will only have to calculate the 
contribution coming from the limits of integration in p° at order 0(1). Indeed, we have: 



dp 



dS p (p°^k) 



dk 



k=ki 



7V 



2 



6 



2klV 



G vp ~ pp (p°,p, x N ) 



Therefore, the expression Hi] gives: 



This result is only for the convolution with the proton-hole spectral function, similarly there 
will be a corresponding term for the neutron-hole spectral function. But the changes are 
minimal. They reduce to change p (protons) by n (neutrons) in the structure functions and 
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5 — > —5. With these changes, we can write the equivalent expression to the above one for 
the convolution with the neutron spectral function as: 

d 3 p M f 



2 M A J d A r J * | J dp»S n h (p",P,k F )G"(p»,p } x N 



+ 



7T 2 5 



7T 



dp { 



k=kt 



G™- pn (p°,p,x N ) H49) 



2Mk F (r)V~ n ^ 11 ' 1 ■ 2k 2 F V dk 

where now G un ~ un (p° , p, xat) is the same as Eq. (144 p but with the replacement of p — >■ n in 
the structure functions. 

Now, we sum Eqs. (1480 and ( )49l) . keeping in mind that now the spectral functions are 
identical as they are evaluated at k = hp. And the result is: 



2 M A I d 3 r 



d 3 p M 
(27r)3£(p) 



dp S h (p°, p, k F ) (G^ p (p , p) + G vn ~ v \p\ p)) 



+ 



+ 



7T 



7T 

2~k\V 



ShOi, P , M p) - G^(/i, p)) 

(G^-^Cp , p) - G up ~ 9p (p°, p)) 



2Mk F (r) V 
2 5 r ^dSnif^k) 



dk 



k=ki 



(50) 



The sum of the structure functions G un ~ un (p°, p) and G up ~ up (jP, p) in the first line of the 
above expression is zero in the limit of considering only light quarks (up- and down-quarks 
approximation), and 



G vn - Vn {p\ p) - G vp ~ Vp {p\ p) = {u v - d v ) 



vp-vpf 



2 4p 2 



M M(p-q) 

where u v = u — u and d v = d — d are the up and down valence distributions, respectively. 



(51) 



Furthermore, in the denominator of the PW relation we have the difference F% A — F% A , 
which is written as: 



f(P* + 2 (p ) 2 -pi) - 2(q ) 2 ((p ) 2 +p 2 ) + 4AW(9°) 2 -9 2 ) • 
F 2 A = = g° — 



2M 2 (g 2 - (g ) 2 ) 



^(p 2 + 2(p ) 2 - p 2 ) - 2(g ) 2 ((p ) 2 + p 2 ) + 4pYp z VW 3 ? )•£,.„ 



2M 2 (g 2 - (q ) 2 ) 

Now, let us define the kinematic factor which goes in front of F 2 like: 



G(p°,p) = 



q 2 (p 2 + 2(p ) 2 - p 2 z ) - 2(g ) 2 ((p ) 2 + p 2 ) + 4p°q°p z ^(q^ 2 - q 2 



2M (q 2 - (q ) 2 ) ■ (p ■ q) 



(52) 



(53) 
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This expression for F 2 is valid for neutrino as well as for antineutrino, therefore, when 
performing the subtraction, we get the following 



Defining 



F% A - F» A = (G(p°,p) (F? - F?)) sl + (G(p°,p) (FT ~ *T)> 5 . 



GT VV {V\V) = G(p°,p) (F? - Ff ) 
GT- pn (p ,P) = G(p°,p) {F™ - F™) 



(54) 



(55) 
(56) 



and using the same procedure as we did for F 1 , one gets 

{c v r 9p {p\p)) t 



'si 



dp S p h (p°,p,k F )G7- up (p",P) 



vp—vpi 



-Sl(fi,p,k F )G2 P up (ii,p) 



5 7T 2 

2VM~k F 

J_t^ r dp0 dS{{p\p,k) 



2Vki 



dk 



g7- Pp (p°,p) + o. 

k=k F \ V 



(57) 



h F J -oo 

where now, the symbols (■ ■ ■) stand to indicate the integrals in d 3 p and d 3 r. And there is 
a very similar formula for the convolution with the hole spectral function for neutrons with 
the known changes (p — >■ n and 5 — > —5). We readily obtain: 



(GT- Dn {p\p)) sn 



+ 



+ 



dp o s n h (p°,p,k F )Gr^ n (p°,p) 

2 -S n h (v,P,k F )G» n -* n (v,p) 



2VMk F 

2 rfpp dSJt{jP,p,k) 



5 7T 



fci 



oo 



k=k F \ V 



(58) 



Therefore, when summing over both the contributions to obtain F% A — F£ A , we get two 
kinds of contributions: the first order one and the second order one (proportional to £). 



F% A - Fl A = (G» p - pp (p°,P)) sfi + (GT- pn (p ,P)) s , 

dp S h (p°,p, k F ) (G?- PP (P°,P) + P "(P°, P)) 

() ' S h (^p,k F ) (GT~ Pn (p ,p) - G»r Pp (p ,p)) 



+ 



+ 



2V Mk F 

5 [» Q dS h (p°,p,k) 
T / (/ /' 



2Vk 



F J -oo 



dk 



{GT- pn ( P °,P) - G7- PP (P°,P)) +o(^- 

k=k F \ V 
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And if one writes now the sum and the difference of functions G2, one obtains (in the up 
and down quarks approximation): 



G^ p (p ,p) + G» n -* n (p°,p) = G(p°,p) (F? - F%> + F™ - F™) = 

G» 2 n -™(p°,p) - G" 2 p - Pp (p°,p) = G(p°,p) {-F»v + F? + FT - n n ) = G(p°,p) x 4x N (u v - d v ) (60) 

Therefore, in this limit, we obtain a dominant non-isoscalarity correction (because the 
isoscalar term is zero) proportional to y. This is exactly the same effect we got when 
calculating F» A - F[ A . 

For the F 3 structure functions for the CC case we have: 

+2A vJ dr J wfm)-J p Sh{p ^ pMn) ip- q ) ( s + 3 } 

= (G7 +PP (P°,P)) SP + (Gr +Pn (p°,p)) sn (61) 

h h 

where in the last step we have defined: 

G u 3 p+Pp (p°, p) = 2A^- P ° q * - q ° Pz (F" p + F 3 Pp ) (62) 
Qz (P-Q) 

G» n+Pn (p°,p) =2A^- P Q j^~Q Pz (F» n + F» n ) (63) 
And now we perform the expansion of the spectral functions around the Fermi momentum: 

{GT 9P {P°,P)) SI = (f_ dp Sl(p ,p,k F )Gr +Pp (p ,& 

} t2 s p h (v,p,k F )Gr pp (^P) 



2V Mk F 



The changes to obtain the G3 convoluted structure function for neutrons are the replacements 
p — >■ n and 5 — > —5, which results 
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Summing the above two equations for the proton and neutron, we get the following: 

F£ A + F* A = I J* dp S h (p°,p, k F ) (Gf + *V,Pl + Gr +Pn (p°,P))\ 



] " 2 -S h ( P ,p, k F ) {GT +9n (^P) - G- p+9p (^P)) 



2V Mkp 

+ (ww F I^tL, (or**V.» - or«W,») + o (A) 2 ) (06, 

If now we write explicitly the combinations of G 3 functions in terms of PDFs considering 
only light quarks (up and down), we get: 

G» p+Pp (p°,p) + G^ n {p\p) = 4G 3 (p ,p) K + d v ) (67) 

where 

G 3 (p°,p)=2A q - P ° q ;- q ° Pz (68) 
Qz (p • q) 

And for the difference we get exactly zero i.e. 

GT +Pn iP°,P) - G 3 W V, p) = (69) 

Therefore, for F 3 , the dominant contribution is the isoscalar one (the one which is not 
proportional to y) and the non-isoscalar contribution (proportional to ^) is zero at first 
order (O (£)). Of course there will be corrections at higher orders, but we are retaining 
only the first order non-isoscalar corrections. 

For neutral currents, in the numerator of the PW relation, we would only have the 
structure function F® A , which can be written as: 

rQA f „ ^ _o„ 1° [ja„ f d3 P M r,JcP/„o^ 



+ 2A— dr / dp S h {p ,p,k F , n ) F 3 U {x N ) 70 

q z J J (2tt) 3 E(p) J_ 00 p ■ q 

^(G^(p^p)) sPh + (G 3 n (p°,p}) s , (71) 



where for the neutral current case, we define 

G^{p\p) = G,{p\p)F^{x N ) (72) 
G^{p\p) = G,{p\p)Fl n {x N ). (73) 

G^(jP,p) is defined in Eq. (|68|) . 
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The expression for F 3 A around the isoscalar condition is: 

F° A = dp°S h (p°,p,k F ) {Gl p (p°,p) + Gl n (p°,p))^ 

+ (wm; Shi ^ kF) ( g T(^p)-g° 3 p (^p))^ 

For the neutral current in the parton model 

if = 2 [(e 2 L (u) - e 2 R (u)) («-«)+ (4(d) - e\{d)) (d - d)} = 2 [«£(„) («-«) + (d) (d - d)] 

(75) 

where 

gl{u) = el{u)-el{u) = \- 2 -si^e w (76) 

^ 2 (d) = e i(d)-4(d) = i-isin 2 % (77) 

and invoking isospin symmetry, we also have: 

Fl n = 2 [«£(«) (d - d) + <?!(d) (« - «)] (78) 

With these, we obtain the expressions for the sum and difference of neutral currents G3- 
functions for protons and neutrons i.e. 

G°V,P) + G° 3 n (p°,p) = 2G 3 (p°,p) {g 2 _(u) + gl(d)) (u v + d v ) (79) 

G° 3 n (p°,p) - G° 3 p (p°,p) = 2 G 3 (p°,p) (g 2 _(d) - gi(u)) (u v - d v ) (80) 

D. Nuclear corrections to PW ratio 

Now using the expressions obtained in the previous section, we may write the ratio of the 
differential scattering cross sections as 

ao NC aa NC 

p— dx dy dx dy 

~ d^ A 

dx dy dx dy 

2xy (1 - f ) [2 [gl{u) + gl{d)) ( dp S h G' 3 (u v (x N ) + d v (x N ))) + O(e 3 A ) 
~ y 2 x 0{e\ A ~* A ) + {l-y~^} 0{et~ VA ) + xy (l - f ) (4 dp* S h G' 3 (u v + d v )) 

Axy (1 - I) (g 2 _(u) + gl{d)) (j^ dp S h G' 3 (u v (x N ) + d v (x N ))) + 0'(ef ) 

= 7 " \ " (81) 

Axy (1 - f) ( dp* S h G' 3 (u v (x N ) + d v (x N ))) + 0'{f A ^ A ) 
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where 



G 



A 



0\el A ) = 2xy (l - f ) 0( e ° 
Q\t vA -~ vA ) = y 2 x 0(ef- pA ) + 



(82) 
(83) 

(84) 



And e indicates that the contribution is of the order O (^). 

Expanding the denominator of Eq. (IHTj) using Taylor series and keeping only the first 



order term in y, we may write the above equation as: 



R- 



Axy (l - |) (g 2 _(u) + g 2 _(d)) /J" dp S h G 3 (u v + d v )\ + 0\e 



i ( r vA—vA\ 



O'ie 



Axy (1 - |) (J^ dp* Sh G' 3 (u v + d v )) 4xy (i _ |) (j^ dp o Sh G> 3 (u v + d v 
g 2 _{u) + g 2 _{d) + 



0'{el A ) 



Axy (1 - |) ( dp° S h G 3 (u v + d v ] 



{g 2 _(u)+g 2 _(d))0'( 



-vA-vA\ 



Axy (1 - §) ( dpo 5, G3 (u v + d v ) 

1 

2 



^ - sin 2 6>w + 5R~ 



15) 



where 5R are the corrections due to the non-isoscalarity of the target and is written as 



5R- = 5R± + 5R~ 2 



(86) 



The first correction 5R l is given by: 



0'{el A ) 



O(e 0A ) 



Axy (1 - I) (J^ dp S h G 3 (u v + d v )) 2 (J^ dp S h G 3 (u v + d v 



2V Mkp 



i-s h ^p,k F )(Gr^Pi-G' 3 op ^,m + ^j^d P od -§ , , (Gr(p°,P)-G' 3 op ( P °,p)) 



k=ki 



2(J^ oo dp S h G 3 (u v + d v : 



(g 2 _(d) - g 2 _(u)) ( &jfc(S h G 3 (u v - d v )) 



n 2V ki, J— 00 r ok 

p J =H F 



k=kp 



G 3 ( y u v d V/ 



J-oc dp S h G 3 (u v + d V/ 



(87) 
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and the second correction 5R 2 is given by: 

,< e vA-vA\ 



(g 2 _(u) + g 2 _{d)) 



- {g 2 _(u) + g 2 _(d)) 



Axy (1 - f ) (J^ dp° S h G' 3 (u v + d v ) 
y 2 x0(e» A -* A ) + {l-y-M^}0(, 



.vA—vA\ 



where 



5R 2A = - (g 2 _(u) + g 2 _(d)) 
5R 2 , 2 = -(g 2 _(u) + g 2 _(d)) 



Axy (1 - |) (J^ dp* S h G' 3 (u v + d v ] 

y 2 xO{el A -* A ) 



5R 2 , + 5R, 2 (81 



Axy (1-|) ( df S h G' 3 (u v + d v ) 



4x2/(1-1) (f^diflSHG^Uv + dJ 



Making the substitution of 0(e\ A ~ DA ) in Eq. ((89]), we get: 



5R- 1 = -(g 2 _(u) + g 2 _(d)) 



y 2M {ww^S h (/,,p, k F ) (G™ - G7-* p ) p0 __ 



- (g 2 _(u) + g 2 _(d)) 
In the parton model, we have: 



With this, we can write: 



M 



«£i = -(0-H+0-(<O) 



_s_ 

2V 



jvffcj [ShGi(u v - d v )) 



k„ J— oo * ok 



(89) 
(90) 



(91) 



[u v -d v )G x {p°,$) (92) 



k=kf 



Gi (u v - d v ) 



2-y 



J^dpOShG'^Uv + dv) 



(93) 



Similarly using G v 2 n vn — G v 2 " p = G2(p°,p) ■ Axn(u v — d v ), we may write 



5R^ 2 = -{g 2 _(u) + g 2 _(d)) 



1 _ _ May 

xy (1 - I) 



X 



2V 



^(S h G 2 x N (u v -d v )) n +^F OQ dp t 



k\ J—oo Qk 



k=kp 



G 2 xn{u v — d v ) 



J^dpVShG'^Uv + d, 



(94) 
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We must point out that when p° = //, Sh((J>,p, kp) is zero [38|, |40| and the imaginary- 
part of the self-energy of the nucleon is also zero, and this has been used in the definition 
(numerator) of the hole spectral function while performing the numerical calculations. 
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FIG. 1. ^-^^h vs y at different x for charged current v^{y^){E Vii = 65 GeV) induced reaction 



in 56 Fe. Dotted line is the base results and for numerical calculations CTEQ 



46|] PDFs at LO 



have been used. Dashed line is the fu 
experimental points are from CDHSW 



1 model at LO. Solid line is full calculation at NLO. The 



3 



(solid circle) and NuTeV Q (sol id square) experiments. 



III. RESULTS AND DISCUSSION 

In this section we present and discuss the results of our numerical calculations. For 
performing the numerical calculations we have used the expressions for the weak structure 
functions and F 3 A for an isoscalar target [lpj ]: 

= 4 / Sr I ^LJL jy S >(P°,P,„W)^ (l + ^) F»(x K ,Q*)(95) 
fjW) = 4/*/|^/Va<*p,*)>^^,SW> (96) 
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FIG. 2. -j; vs y at different x for neutral current u fjL (u lJL )(E Uii = 65 GeV) induced reaction in 
56 Fe. Lines have the same meaning as in figHJ 



and for a nonisoscalar target, the expressions for F A (x) and F A (x) are given by llj : 



F 2 A (x A ,Q 2 



d 3 r 



d 3 p M 
(2tt)3£( p ) 



+ / dp S^ utron (p ,p,k F;n )F^ utr<m (x N ,Q 2 ) 



dp° sr ton ( P °, P , k F , p )Fr ton (x N , q 2 ) 

2x N p 2 x 



x 



l + w 

Xn V Mun 



(97) 



F 3 A (x A ,Q 2 



d 3 r 



d?p M 
(27r)3£(p) 



Hp 



dp S Proton {p0 ^ ^ k ^jf"*"^ (?) 



dp S™ utron (p°, p, k F<n )F£ eutron (x N , Q' 2 ) 



P°7 ~ Pz 

(p° - Vzi)i 



where 



7 = ^ = l 1+ Q 2 



4M 2 X 2X 1/2 



Xm 



Q 2 



2(p°q° -p z q z 



(98) 



(99) 



Here Fj'g * ' 1 and F^ utron are the dimensionless structure functions for the free proton and 



2,3 



the free neutron respectively. These structure functions are calculated with target mass 
correction (TMC) [43| and CTEQ6.6 parton distribution functions (PDFs) at the Leading- 
Order (LO) [461 ] . Fermi motion and nucleon binding are implemented through the use of a 
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nucleoli spectral function. This is our base result. We then also include pion and rho cloud 
contributions in F^ following the model of Ref. 40J and shadowing corrections in F^ and 



441 ]. which is our full calculation (Total). Using them we have obtained the charged 



current differential scattering cross sections, the expression for which is given by Eq. (|25p . 
For the neutral current the expression would remain the same with the changes given in 
Eq. (|26|) . Using these cross sections, we evaluate the ratio 



R~ 



rf 2 o-(i/ M N-^v^ X) 
dxdy 



d 2 a(u,, N^tu^ X) 
dxdy 



cPaju^ N^-fi' X) 
dxdy 



dZatyVy, AT->/x+ X) 
dxdy 



(100) 



Results for the charged current induced (anti)neutrino-iron differential cross section have 
been discussed in detail in Ref. [loj], where the comparisons have been made with the exper- 
imental data of NuTeV and CDHSW which are corrected for isoscalar iron target, however, 
for the completeness we are showing here the results for the charged and neutral current 
induced (anti)neutrino-iron differential cross section in FigsJH and [2] respectively. In the 
case of neutral current we have performed numerical calculations at LO. We observe that 
medium effects are important in the study of differential scattering cross sections. For ex- 
ample, when the calculations are performed for charged current neutrino induced process at 
LO, medium effects like Pauli blocking, Fermi motion and nucleon correlations which are 
taken into account for our base calculations, results in the reduction of the differential cross 



E=80 OeV 



0.25 — s y 



x=0. 1 25 
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x=0.55 

Isoscalar 



FIG. 3. Paschos and Wolfenstein ratio R 



d 2 o[v tl N-Ufy, X) 
dxdy 



d^a^fj, N-^Vfj, X) 
dxdy 



in 56 Fe for non- 



d 2 <r(y (J N-*ii- X) d 2 <r(ii M JV->p+ X) 

dxdy dxdy 

isoscalar case. R~ is calculated for (anti)neutrino energy(E) of 80GeV and at different Bjorken x. 



The solid line is the result when one treats 56 Fe to be an isoscalar target. 
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FIG. 4. Nonisoscalar correction(<5i? ) vs y at different values of x for (anti) neutrino energy(E) of 
80GeV. 



E=80 GeV 



■ ■ ■ ■ x=0.2 

■ x=0.3 
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- ■ ■ — x=0.5 
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FIG. 5. sin 2 8 w vs y at different values of x in 56 Fe treating it to be non-isoscalar nuclear target 
for (anti)neutrino energy(E) of 80GeV. 

section which changes by 3 — 4% at low y and increases to 10 — 11% at high value of y for the 
studied region of x from the free case. And when pion and rho cloud contributions as well 
as shadowing effect are also taken into account there is a further change of about 9% at low 
y and this change decreases with the increase in y for x=0.2-0.3. This difference becomes 
smaller with the increase in the value of x, for example at x=0.5-0.6, this becomes less than 
1% for all values of y. Similarly, in the case of antineutrino induced charged current process 
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FIG. 6. sin 2 6 w vs x in 56 Fe treating it to be non-isoscalar nuclear target. The results are shown 
at different values of Q 2 and (anti) neutrino energy E. 

the change in the base results from the free nucleon scattering process is around 10 — 12% 
at x=0.2-0.3 for studied region of y, which gets further modulated by 12% at low values of 
y and significantly increases with the increase in the value of y when pion and rho cloud 
contributions and shadowing effects are incorporated. At higher values of x, the difference 
in the results obtained using the base calculation and the full calculation is negligibly small 
and the difference between the base results and free case is same as in the case of neutrino. 
Furthermore, calculations performed at NLO make the results better. In the case of neutral 
current we have performed calculations at LO and observe that the difference between the 
base results and the results for the free case (not shown here) is 10 — 12% at x=0.2-0.3 for 
neutrino as well as antineutrino induced processes. When pion and rho cloud contributions 
and shadowing effects are taken into account results from base changes to 17 — 25% for 
neutrino case and for antineutrino induced process the change is 30 — 40% at x=0.2-0.3 for 
all values of y. At higher x, when Fermi motion, Pauli blocking and nucleon correlation 
have been taken into account, the results changes by 4 — 10% from the free results at all 
the values of y for both neutrino as well as antineutrino induced processes. When pion and 
rho cloud contributions as well as shadowing effects are also considered, there is a further 
reduction of 2 — 3% for neutrino induced process and 4 — 9% for antineutrino case for all 
values of y. Therefore, we observe that in the study of charged current and neutral current 
differential scattering cross sections medium effects are important. 
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The PW ratio(i?~), given by Eq. fllOOp . has been shown in Fig J3] using the numerical 
values from charged and neutral current differential scattering cross sections in iron, treating 
it to be isoscalar as well as nonisoscalar nuclear targets, at different values of x for a fixed 
value of neutrino/antineutrino energy E u ^=80 GeV. We have incorporated Fermi motion, 
Pauli blocking and nucleon correlation while calculating R~. We find that R~ is almost 
independent of x and y for an isoscalar target, while for the nonisoscalar target there is x 
as well as y dependence. We find that the nonisoscalarity and medium effect increases with 
the increase in the value of x, for example in the mid region of x it is about 14 — 15% for 
low values of y and 3 — 4% at higher values of y. Therefore, for nonisoscalar targets like 
iron, medium effects and excess of neutrons strongly affect R~ and hence sin 2 9w While 
for isoscalar target like carbon(not shown here) medium effects cancels out and extracted 
value of sin 2 9w from PW ratio is in complete agreement with the global value. When the 
contributions of pion and rho clouds are also taken into account, we find that these changes 
do not bring any difference in the results of PW ratio as obtained with our base calculations. 
Since shadowing is a low x phenomena therefore we have not considered it in our present 
study of PW ratio [S|. To see the effect of nonisoscalarity in the iron target, we have plotted 
SR~ using Eq. fl86|) at E v (p\ = 80GeV in FigJH We find that the effect of nonisoscalarity is 
large at low y and high values of x which decreases with the increase in the value of y as well 
as this effect is smaller at low values of x. Hence there is a nonisoscalarity dependence in the 
determination of sin 2 9w- Using the results of R~ from Eq. ( 11001) and SR~ from Eq. ( l86l . 
we have obtained sin 2 9w using Eq. ( 1851) and presented the results in Fig. El The above 
calculated value of weak mixing angle is now corrected for isoscalar target where we should 
also keep in mind that the medium effects are still present. We find that due to medium 
effects sin 2 9w is different from the global fit, and this difference is ~ 7% when evaluated 
for low value of y at x=0.2 and this decreases to 1% at high values of y, while this change is 
~ 9% when calculated for low y at x=0.6 and this reduces to 2% at high values of y. Thus 
we observe that nonisoscalarity as well as medium effects like Pauli blocking, binding energy 
and Fermi motion are important while extracting sin 2 9\y- To see the effect of neutrino 
energy E v and Q 2 dependence on sin 2 9w, we have plotted in Fig® sin 2 9w as a function 
of x at various values of E v and Q 2 and we observe that at £^=80 GeV and Q 2 =25 GeV 2 
it is almost close to the standard value at low value of x and the value of sin 2 9w changes 
significantly with E u / P ), Q 2 and x. Therefore, while looking at the NuTeV results it is also 
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important at what values of x, E u and Q 2 , the analysis were performed due to the fact that 
a wide range of these variables could change considerably the value of sin 2 9w- 



IV. CONCLUSIONS 

To summarize our results, in this work we have studied the effects of the nuclear-medium 
on the structure functions F^(x,Q 2 ) and F^(x,Q 2 ) in the iron nucleus treating it to be 
isoscalar nuclear target and then we have made non-isoscalarity corrections. We have used 
many-body theory to describe the spectral function of a nucleon in an infinite nuclear- 
medium. The spectral function takes into account the Fermi motion and the binding energy 
of the nucleons. The spectral function also includes nucleon correlations. Then to apply it to 
the case of finite nucleus we have used local- density approximation. Target mass corrections 
have also been considered. We have used CTEQ 46] PDFs in the numerical evaluation. We 
have taken into consideration the effects of mesonic degrees of freedom, of shadowing and 
antishadowing for the calculation of F% and for shadowing and antishadowing effects 
only. 

With these structure functions, we have evaluated numerically the differential scattering 
cross sections for the charged and neutral current induced (anti)neutrino interactions on the 
iron target. These differential scattering cross sections are then used to study the nuclear 
medium effects and nonisoscalarity correction in the extraction of weak mixing angle sin 2 ^ 
using Paschos-Wolfenstein(PW) relation. 

Beginning from a model for non-symmetric nuclear matter, we have expanded the nu- 
clear hadron tensor in an isoscalar part plus nuclear corrections which are, roughly speak- 
ing, proportional to the difference between the neutron and the proton density profiles. 
We have performed this expansion for all the nuclear structure functions which appear in 
the Paschos-Wolfenstein ratio for differential cross sections and obtained the pure isoscalar 
Paschos-Wolfenstein result plus nuclear corrections which depend on the phase space Bjorken 
variables (x,y). 

We have also observed that the Paschos-Wolfenstein ratio for isoscalar nuclear structure 
functions is almost independent of (x,y), as it should be. However, when one considers the 
model for non-symmetric nuclei, such as 56 Fe, one obtains a PW ratio which depends on the 
phase space variables. Furthermore, we have evaluated the first order nuclear correction to 
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the isoscalar PW ratio and found it to be important for a wide range of the phase space. 
If this nuclear correction can explain the NuTeV anomaly is a question that should be 
answered if one could weight our results for 8R~ by the relative amount of events in every 
(x, y) bin, under the assumption of assigning all the correction to sin 2 6w 

Another important point that should be taken into account is the fact that we have worked 
in the up and down quarks approximation, therefore, we have neglected the contribution 
coming from heavier flavors. We leave this subject for future studies. 
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